Adaptive time-stepping for incompressible flow part I: scalar advection-diffusion

Even the simplest advection-diffusion problems can exhibit multiple time scales. This means that robust variable step time integrators are a prerequisite if such problems are to be efficiently solved computationally. The performance of the second order trapezoid rule using an explicit Adams–Bashforth method for error control is assessed in this work. This combination is particularly well suited to long time integration of advection-dominated problems. Herein it is shown that a stabilized implementation of the trapezoid rule leads to a very effective integrator in other situations: specifically diffusion problems with rough initial data; and general advection-diffusion problems with different physical time scales governing the system evolution

Adaptive time-stepping for incompressible flow. Part II: Navier-Stokes equations

We outline a new class of robust and efficient methods for solving the Navier- Stokes equations. We describe a general solution strategy that has two basic building blocks: an implicit time integrator using a stabilized trapezoid rule with an explicit Adams-Bashforth method for error control, and a robust Krylov subspace solver for the spatially discretized system. We present numerical experiments illustrating the potential of our approach. © 2010 Society for Industrial and Applied Mathematics

Robust a posteriori error estimators for mixed approximation of nearly incompressible elasticity

This paper is concerned with the analysis and implementation of robust finite element approximation methods for mixed formulations of linear elasticity problems where the elastic solid is almost incompressible. Several novel a posteriori error estimators for the energy norm of the finite element error are proposed and analysed. We establish upper and lower bounds for the energy error in terms of the proposed error estimators and prove that the constants in the bounds are independent of the Lam\'{e} coefficients: thus the proposed estimators are robust in the incompressible limit. Numerical results are presented that validate the theoretical estimates. The software used to generate these results is available online.Comment: 23 pages, 9 figure

A preconditioner for the 3D Oseen equations

We describe a preconditioner for the linearised incompressible Navier-Stokes equations (the Oseen equations) which requires as components only a preconditioner/solver for each of a discrete Laplacian and a discrete advection-diffusion operator. With this preconditioner, convergence of an iterative method such as GMRES is independent of the mesh size and depends only mildly on the viscosity parameter (the inverse Reynolds number). Thus when the component preconditioner/solvers are effective on their respective subproblems (as one expects with an appropriate multigrid cycle for instance) a fast Oseen solver results

Iterative Methods for Problems in Computational Fluid Dynamics

We discuss iterative methods for solving the algebraic systems of equations arising from linearization and discretization of primitive variable formulations of the incompressible Navier-Stokes equations. Implicit discretization in time leads to a coupled but linear system of partial differential equations at each time step, and discretization in space then produces a series of linear algebraic systems. We give an overview of commonly used time and space discretization techniques, and we discuss a variety of algorithmic strategies for solving the resulting systems of equations. The emphasis is on preconditioning techniques, which can be combined with Krylov subspace iterative methods. In many cases the solution of subsidiary problems such as the discrete convection-diffusion equation and the discrete Stokes equations plays a crucial role. We examine iterative techniques for these problems and show how they can be integrated into effective solution algorithms for the Navier-Stokes equations